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I. Introduction 

In this note we give a new and, we believe, simpler proof of a recent result of Laptev 
and Weidl. It is concerned with Lieb-Thirring inequalities for matrix valued Schrodinger 
operators of the type 

H = -^®I + V{x) (1) 

acting on L^(]R;C^). The potential V{x) is a negative definite hermitean N x N matrix. 
We assume that its matrix elements are smooth functions of compact support, say in the 
interval [—a, a]. The operator H has finitely many negative eigenvalues, which, counting 
multiplicities, we denote by — Aj j = 1, . . . , L. 
The following theorem was proved in [LW]. 

Theorem 1 

With the above assumptions on V the following inequality holds 

EAf <-|/TV(y(xf)dx. (2) 



From Weyl's law on the distribution of eigenvalues it is seen that this inequality is 
best possible. For the case where the potential is a scalar function this result was already 
proved in [LT] where it was realized that (2) follows from trace identities. 

The matrix case, however, is important, since inequality (2) of Laptev and Weidl is 
the starting point for deriving sharp Lieb-Thirring inequalities in higher dimensions. In 
particular, the argument of [AL] applies also in this case and yields sharp Lieb-Thirring 
inequalities for the sum of powers of eigenvalues where the power is larger than 3/2. For the 
details we refer the reader to the original paper [LW] where a collection of beautiful results 
is presented. Their proof of Theorem 1 which corresponds to formula (2.1) in their paper 
is patterned after the proof of [BF] (see also [FZ]) and is fairly involved. Inequality (2) is 
derived from a trace identity, which in turn is a special case of a whole family of identities 
that express conservation laws of the Korteweg-de Vries equation. The derivation of these 
trace identities uses nontrivial results about scattering theory on the line and Laptev and 
Weidl prove these afresh for the matrix case. Since (2) is the central result in [LW] and of 
independent interest, it is of value to have a different, more elementary and more direct 
proof. It relies on the 'commutation method' and some elementary facts from the calculus 
of variations. 



The 'commutation method' has a fairly long history, some versions of it were already 
known to Darboux [DG] and Jacobi[J]. Its modern appearance seems to be due to Crum 
[C] . For a rigorous discussion of these issues we refer to the papers of [G] and [DP] . In the 
latter more examples of the usefulness of this method are presented. Another work, closer 
to the spirit of ours, is the one of Schmincke [S] who uses the commutation method to 
prove that 

5^Af^>-^|y(x)dx 

for scalar potentials. This result was extended in [LW] to the matrix case which can also 
be obtained using the methods of the present work. This inequality should be contrasted 
with 

j2^y"<-llv{x)dx 

obtained in [HLT] for the scalar case and in [HLW] for the matrix case. Both inequalities 
are sharp in the sense that the constants cannot be improved. 

To illustrate the ideas we give a short proof of Theorem 1 for the case where V^ is a 
scalar potential, thereby recovering the result in [LT] . This sets the stage for the proof of 
the matrix case in the following section. While it is certainly possible to prove Theorem 1 
under fairly general conditions on the potential, we refrain from doing so. It would clutter 
the simple argument with technical details. 



II. The scalar case 

Let — Ai be the lowest eigenvalue of the Schrodinger operator (1) of Section I with 
a scalar potential. It is well known that this eigenvalue is not degenerate and the cor- 
responding eigenfunction (pi can be chosen to be strictly positive. Moreover, outside the 
range of the potential we have 



const. e vAi^a;^ if x > a, 
const. e^'^^, if X < —a. 



^^(^^)^^consl.e^^^ , '}^^>^^ (1) 



Thus the function 

Fi^) = ^ , (2) 

(pi{x) 

is defined and satisfies the Riccati equation 

F' + F^ = V + Xi , (3) 
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together with the conditions 

F(x) = (-/^' ^J^>"' (4) 

I vAi, if a; < -a. 

A simple computation shows that the Hamiltonian H can be written as 

H = D*D-Xi, (5) 

where 

^-^-f. (6) 

and 

D* = -^-F. (7) 

da: 

It is a general fact [DP][G] that the operators D*D and DD* on L^(]R) have the same 

spectrum with the possible exception of the zero eigenvalue. Note that D*D has a zero 

eigenvalue which corresponds to the ground state of H. The operator DD* does not have a 

zero eigenvalue. This follows from the fact that the corresponding eigenfunction t(j satisfies 

i;' = -Fi; , (8) 

and hence '^(x) = const. /(^i(x) which grows exponentially and is not normalizable. Thus 
the new Schrodinger operator 

d^ d^ 

H = DD* -Xi = -—r-F' + F^-Xi = -—r + V-2F'. (9) 

dx^ dx^ 

has, except for the eigenvalue — Ai, precisely the same eigenvalues as H. Also note that 
the potential V — 2F' is smooth and has support in the same interval as the potential V. 
Next, we compute using the Riccati equation (3) 

f{V - 2F'fdx = f V^dx + 4 /"(Ai - F^)F'dx . 

The last term can be computed explicitly using (4) and we obtain 

f{V - 2F'fdx = f V^dx - y Af' . (10) 



Thus, 



t >^f - T(, I "'■^ = t ^f - 1 /<^ - 2^''''^^- ■ '"' 



fe=l -^ k=2 
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and the Schrodinger operator with the potential V — 2F' has precisely the eigenvalues 
— A2, . . . , — A^. Continuing this process we remove one eigenvalue after another. After the 
last one is removed a manifestly negative quantity is left over, and this proves Theorem 1 
in the scalar case. I 



III. The matrix case 

The proof of Theorem 1 is patterned after the scalar case. In addition to the usual 
eigenvalue equation for H in (1) of Section I 

-(l)"{x) + V{x)(i){x) = -X(i){x) (1) 

we consider the following matrix version for a,n N x N matrix M(x), 

-M"{x) + V{x)M{x) = -XM{x) (2) 

The following Lemma is central. 

Lemma 2 

Assume that —A is the ground state energy of H and let (j) he any solution of the 
differential equation (1) with 

(l){x) = e^^^M for X < —a 

where j^ u E C^ is constant. In particular, we do not require that (p is normalizable. Then 
(l){x) never vanishes. Moreover, the ground state energy is at most N-fold degenerate. 

Proof: Suppose there exists a point xq with (I){xq) = 0. Consider the continuous 

function 

)(x), if X < xo 



"0, n X > xq. 

Clearly, this function does not vanish identically and is square integrable. A simple inte- 
gration by parts calculation shows that 
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and thus cj) is a, ground state and must be a solution of the Schrodinger equation (1) which 
is an ordinary differential equation of second order. Here ( , ) denotes the inner product 
on L^(]R, C'^). Since cp vanishes to the right of xq the solution must vanish everywhere, 
which is a contradiction. The last statement of the lemma is an immediate consequence 
of this. ■ 

Remark: The above Lemma clearly generalizes to potentials that do not have compact 
support but decay, e.g., exponentially, at infinity. 

Consider any matrix solution M(x) of the differential equation (2) subject to the 
condition 

M{x) = e^^A for x < -a , (3) 

where A is a nonsingular matrix. By the previous Lemma 1, any solution of (1) that decays 
exponentially must be a linear combination of the column vectors of M{x). In particular, 
the ground states themselves must be linear combinations of the column vectors of M{x). 
Also by Lemma 2 we know that the matrix M(x) must be invertible for every x G M. 
Hence it makes sense to define 

F{x) = M-\x)M'{x) . (4) 

The following Lemma 3 states all we need to know about F{x). The number K below 
denotes the degeneracy of the ground state energy. We have that K < N hy Lemma 2. 

Lemma 3 

The matrix F{x) is hermitean for every x & M., independent of the choice of A and 
satisfies the matrix Riccati equation 

F' + F'^ -V ^XI . (5) 

Moreover, for x < —a 

F{x) = VXI , (6) 

and for x > a, the eigenvectors of F{x) are independent of x and and its eigenvalues 
decay exponentially fast to the K fold eigenvalue —\fX and the N — K fold eigenvalue vA 
respectively. 
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Proof: Consider any two matrix solutions of (2), Mi{x) and M2{x). From the Wron- 
skian identity ( with * denoting adjoint) 

-^ (MiixYM^ix) - Ml'{x)M2{x)) = , (7) 

da; 

one obtains 

(Mi(x)*M^(x) - M*'(x)M2(x)) = const. . (8) 

First we set Mi = M2 = M and assuming the initial condition (3) we get that 

M{xyM\x) = M*'{x)M{x) , (9) 

which yields the hermicity of F. If we set Fi = M[M^ and F2 = M2M2 where Mi and 
M2 satisfy (2) and (3) for possibly two different, nonsingular matrices Ai and A2 we get 
from (8) that Fi = F2. An elementary computation yields (5) and (6). 

Fix xq > a. For x > xq the potential vanishes and the matrix M{x) is given by 

M{x) = cosh(VA(x - xo))M{xq) + —= sinh(v^(x - xq))M' {xq) , (10) 

V A 

and hence F{x) is given by 

F{x) = \/\ ('v^tanh(v^(x - xo))I + F{xo)\ (\fXI + tanh(v^(a; - xo))F{xq)\ 

(11) 
From this it follows that the eigenvectors of F{x) do not depend on x and since 

F{x) exists for all x we must have — vA/ < F{xq) < yXI. It follows from (10) that the 

bound states are precisely those solutions (l){x) of the differential equation (1) that decay 

exponentially in both directions and that are of the form (t){x) = M{x)M~^{xo)u where 

u is an eigenvector of F{xo) with eigenvalue — V^- Thus F{xo) has the K fold eigenvalue 

— vA and all the other eigenvalues Vj satisfy the inequality — vA < Vj < vA. From (11) 

we see — V^ is a K fold degenerate eigenvalue of F{x) for all x > xq and that all the other 

eigenvalues converge exponentially fast to vA. ■ 

Proof of Theorem 1 in the matrix case: From the Riccati equation (5) we get that 

H + XiI = D*D, (13) 

where 

D*=(-^^I-f) , (14) 



and 

Z>=(A«,_f). (15) 

Clearly 

D(j) = (16) 

for any ground state (p. Moreover, 

D*t(j = (17) 

has no nontrivial normalizable solution on M since F = yXI for x < —a. 
Thus the operator 

H' := DD* - XJ (18) 

has precisely the eigenvalues A^+i, • • • , A^- A calculation shows that 

d^ 
H' ^-—®I + Vix)-2F'{x) (19) 

where the potential 

V{x)-2F'{x) (20) 

is smooth and decays exponentially fast at infinity by Lemma 3. One easily computes 
using (5) that 

f Tr{{V- 2F'f) dx 

= f Tr (V^) dx-4 f Tr {{F^ - Ai/) F') dx , (21) 

Jr Jr 

which can be integrated to yield 

-^Tt {F\x)) \tZ + 4AiTV {F{x)) \+Z ■ (22) 



By Lemma 2 this equals 

~ 3 



'4KXr ■ (23) 



Again, we have shown that 

L „ ^ L 



5 Af - 1 /^ T,V^)d. . E .f - I^ /.-•((- - -0^) d. 



and the Schrodinger operator with the potential V — 2F' has precisely the eigenvalues 
Ak+1, • • • , Al- 



The potential V — 2F' decays exponentially but, unfortunately, does not have compact 
support and hence the second step in the scalar case, i.e., the removal of the next eigenvalue, 
cannot be taken directly. However, the following approximation argument can be used to 
circumvent this difficulty. Cutting off the potential V — 2F' sufficiently far out we are 
left with a new potential V — 2F^ which has compact support and whose eigenvalues are 

numbers —hk+i, ■ ■ ■ — A*i which can be made to be as close to the old ones — A/c+i, ■ Al 

as we please. The cutoff might cause some new eigenvalues to appear, but all of those can 
be made to be as close to the continuum, i.e., as close to as we please. Removing the 
ground state eigenvalue of this new potential V — 2F^ yields 



t, Mf -|/Ti-((y-2F3-)dx 



j=K+l 

E ^f-^[Tr{{V-2F^-2Gr)dx. 

j=K+P+l ^ 

Here P denotes the degeneracy of Hk+i and G plays the same role for V — 2F^ as F does for 
V. Although tempting, one cannot remove the cutoff in this formula since the two terms 
on the right side are not separately continuous. E.g., the degeneracy of the eigenvalue 
Hk+1 is not necessarily the same as the degeneracy of the eigenvalue Xk+i- Nevertheless 
we have the following 



t Af -^/Tr((V--2F'Hcte 



j=K+l 

E ^f -^/Tr((y-2F^-2Gr)dx +ei, 

where ei is the error in the eigenvalues and the potential integral due to the cutoff in the 
potential. Again, the new potential has exponential decay. 

By repeating the cutting and removing procedure finitely many, say s < L times, we 
end up with 

E ^T - ^ / Tr(V^')dx = -^ / T^{W^)dx + ei + • ■ ■ + e, , (25) 



where the e^ denotes the error stemming from the cutoff at the i-th step and W is the 
resulting potential. In particular (25) implies that 



t\'"-;|/T^(^')d-£-i + 



which we can make as small as we please. I 
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